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We present a combinatorial method of constructing solutions to the normal ordering 
of boson operators. Generalizations of standard combinatorial notions - the Stirling and 
Bell numbers, Bell polynomials and Dobinski relations - lead to calculational tools which 
allow to find explicitly normally ordered forms for a large class of operator functions. 

Key words: boson normal ordering, coherent states, combinatorics 

1 Preface 

Solutions to mathematical problems which we would call nowadays examples 
of solutions to the normal ordering problem have been known from a long time. 
Transformation of a differential operator (a; 3^)" into a form in which all derivatives 
stand to the right with respect to all powers of x, or the Baker-Campbell-Hausdorff 
formula which allows to disentangle the exponential operator exp (A + B) for oper- 
ators A and B satisfying commutation relations [A, B] — C, [A, C] = 0, [B, C] = 
provide us with well known examples. In the framework of the second quantization 
scheme normal ordering of expressions involving ladder operators, i.e. transfor- 
mation of such expressions into a form in which all creation operators are moved 
to the left with respect to all annihilation operators using appropriate commuta- 
tion relations, plays an important role as it simplifies evaluation of the vacuum 
and coherent state matrix elements in quantum field theory, many body quantum 
physics and quantum optics. A standard way of reaching the goal is provided by 
the Wick theorem - an universal tool, unfortunately ineffective when one wants to 
investigate general structure of solutions for a vast majority of physical problems. 
Thus it is highly desirable to find manageable formulas or at least some guiding 
principles leading to the solutions of general normal ordering problem formulated 
as follows: consider an operator F{a'' , a) depending on canonical boson ladder oper- 
ators [a, a^] = 1 and search for solutions G{a\ a) of a functional operator equation 



F{a\a) = N[F{a\a)] =: G{a\a) 



(1) 
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where normally ordered operator N[F{a\a)] is obtained using nontrivial commu- 
tation relations while under the symbol : : the ladder operators are considered as 
commuting ones. Research on general normal ordering problem started from semi- 
nal paper of Cahill and Glauber U| who considered operators bounded with respect 
to the Hilbert-Schmidt norm, expanded them in terms of ordered products of lad- 
der operators and found integral representations for coefhcients of such expansions. 
The method, however, allowed to solve Q explicitly only for the simplest exam- 
ples. Investigations were pushed forward when Navon [2] and Katriel |3J showed 
that the normal ordering problem is situated firmly in the area of combinatorics 
and that methods of combinatorial analysis do appear effective in finding solutions 
of Since that time many examples have been studied and normally ordered 
forms of operators depending on (at)''a^ characterized by integer powers r and s 
or depending on q{a^)a + v{a^) with arbitrary functions q and v were found 0] - 
[7]. Here we are going to present basic concepts of our approach and to sketch 
combinatorial methods useful in order to solve the boson normal ordering problem. 



2 Stirling and Bell numbers, Dobinski relations - the generic example 

Considering canonical boson ladder operators a and we can ask what is 
normally ordered expression for the n-th power of the number operator N — a'^a. 
Obviously such an expression may be written down as 

n 

iV" = (a^a)" = J2 ^('^' k){a^)''a'' (2) 

but in order to solve the problem one has to determine ^(n, fc)'s. This is easy to 
do using mathematical induction - the numbers S{n, k) satisfy the recurrence 

S{n + l,k) = kS{n,k) + S{n,k-l) (3) 

which defines the Stirling numbers of the second kind - positive integers which 
appear in enumerative combinatorics ^) or in classical problems: 



^-^) ^Y.^S{n,k)x^-^, ^ Z^S{n,k)x^ (4) 



where x- = x ■ {x — 1) ■ ... ■ {x — k + 1) denotes the falling factorial {x- =1). The 
second equation in called the Stirling transform, describes transformation of 
basis in the space of polynomials so the Stirling numbers of the second kind may be 
also interpreted as coefficients leading from {x— to {x"}^^q, while the Stirling 
numbers of the first kind provide us with coefficients of the reverse transformation. 
The basic tool extensively used in our further considerations is the Dobinski for- 
mula. In order to derive it let us apply to the exponential e^. This gives 

^) The Stirling numbers of the second kind count the number of ways of putting n different 
objects into k identical containers none leaving empty 
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oo /c ^ 

fc=0 ■ k=l 

namely the Dobinski formula, which is a remarkable relation between some polyno- 
mials B{n, x) (called Bell, or exponential, polynomials lOj) and products of infinite 
series. List of implications of jSjl reads 

— positive integers Bn — B{n, 1), called Bell numbers, can be written for any n as 
infinite sums of fractions 

— the Stirling numbers of the second kind can be explicitly calculated from 

• j=i \ / 

— a closed formula for exponential generating function of the Bell polynomials can 
be obtained 

oo ,„ oo oo \n 

GiX,x) = ^ B(n,x)^ = e- E E fc"frtr 

n=0 n=Ofc=0 (o\ 

Using the properties of standard coherent states (defined as eigenstates of the anni- 
hilation operator, a\z) = z\z)) we conclude that the diagonal coherent state matrix 
elements of N"^ generate Bell polynomials: 



{z\{a^aY\z) = (z| S{n, k)a^''a''\z) = ^ S{n, k)\z\''^ = B(n, |z|2). (9) 

fe=l k=l 

Moreover, expanding the exponential e^°*° (with A being an arbitrary parameter) 
and taking the diagonal coherent state matrix element we have 

{z\e^^'^\z) = g {z\{aUr\z)^ = g B{n, jzp)^, (10) 

n— ■ n— 

which means that the diagonal coherent state matrix element of e'^°*° yields the 
exponential generating function of the Bell polynomials {z\e^°' °-\z) — e'^' 
Relation between the considerations above and the normal ordering problems comes 
from an important property of the coherent state representation (ll) : 
If for an arbitrary operator F{a,a^) we have 

{z'\F{a,a^)\z)^{z'\z) G{z'\z) (11) 



3 



p. Blasmk et al. 



then the normally ordered form of F(a,a^) is given by 

A/'[F(a,a^)] = : G(a^ a) : . (12) 

Consequently, the knowledge of the generating function of the Bell polynomials 
lead to 



= : e'^^'^^'^'-^) :, (13) 



i.e. the normally ordered form of e^*^*". We shall generalize this result, using the 
methods just described, to a quite large class of operator functions. 

3 Normal Ordering of Homogeneous Boson Polynomials 

Let us treat the normal ordering of expressions in boson ladder operators called 
homogeneous boson polynomials and given as linear combinations of monomials 

{a^'Y^'a'" ...{a''Y'a''ia''Y'a'\ (14) 

all having the same difference between the total number of creation and annihilation 
operators which we call excess, denoted by c? X]f=i(''2 ~ ^i) assumed to be 
nonnegative. It means that we consider operators which normal form is 



N 



with Nq — si, N — ■Sj E^nd Q!fc's to be specified for any initially given linear 

combination of monomials (|14|l |H],|1|. 

Calculating the normally ordered form of (H^)" we get 



nN 



{HiY^ia^T'^ ^ Si{n,k) {a^fa\ (16) 

k=No 

with Sa{n, k) to be determined. These coefficients obviously generalize the Stirling 
numbers, provide us with generalized Bell polynomials and Bell numbers 

nN nN 

Btin,x)^ E Si{n,k)x'', Bi{n) = Bi{n,l) ^ E Si{n,k) (17) 

k=No k=No 

and share general properties of the standard Stirling and Bell numbers |Hl,|ni- 
— satisfy the recurrence relation generalizing ^ 

St{n + l,k)= E"' E p )i^d + k-l+pr-Si{n,k-l+p), (18) 
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determine connection between two sets of polynomials 

n N nN 
1=1 k=No k=No 

lead to generalized Dobinski-type relations 

I N 



(19) 



E 

1=0 



n E ii + i^^-m- 



.i=l k=No 



(20) 



Generalized Stirling numbers and their generating functions are explicitly given by 



k / I. \ n N 

si{n.k)^^Y.[ ■ i-^f-' n E (i + - mK 

- J J i=l l=No 



3=0 

oo k / h. 



n=k 



k{ ^ \ J 



n N 

n E akil + ii-l)dy 

i=l k=No 



A" 



(21) 



from which it is seen that generalized Stirling numbers can be expressed in terms 

of ordinary Stirling numbers S{m, I) with m — I, ■ ■ ■ , n{N — Nq + 1). 

If we introduce exponential generating functions of generalized Bell polynomials 



G^(A,:^) = E^'("'^) — 



(22) 



n=0 



then we can derive relations between coherent state matrix elements and suitable 
generating functions. Analogously to (|10|l we have 



A" 



((z*)^A)" 



(23) 



which yields (z|e^^'|z) = ((z*)'*A, \z\^) and 

° = : ) A,a'a) 



(24) 



Inserting the generalized Dobihski relation H2U|I into the generating function (|22|) 
one obtains the latter as 



oo oo 



G^(A,x)=e--EE 



n=0 1=0 

and, after a change of summation order, as 



n E (i + i^-m- 

1=1 k=No 



A" 



OO / C)0 

1=0 ' n=0 



N 



n E i^ + i^-m- 



,i=l k=No 
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(25) 



(26) 
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The series H26(l is convergent either if A'o = = 1 which means that we are 
treating operators including arbitrary number of creators and only one annihilator 
or if d = 0, which corresponds to depending on the number operator. The 
first case may be illustrated on the example of {a'^Ya, and generalized to the 
normal ordering problems of operators involving expressions being linear in a and 
depending on arbitrary functions (/(a^) and v{a^) 

-:g(A,at)e[^(^'"')-°1'^ : , (27) 

where, according to the substitution theorem, both functions g, T are unique solu- 
tions of linear differential equations determined by q and v. Coherent state matrix 
elements of (|27|) give generating functions of Shcffer polynomials . For the second 
case the internal sum becomes exponential of a polynomial in I and the external 
sum converges if aj^X < 0. This leads, [7], to 

oo / ^ \ i 

G^(A,x) = e--^exp A ^Jf (28) 

1=0 \ k=No / 

Combinatorial generating functions (|22|l . introduced as formal power series, may 
be divergent but if we require their interpretation as physically meaningful matrix 
elements then we have to give them analytical meaning, at least in the sense of 
asymptotic expansions. It is our conjecture that in such a case the Pade method 
of generalized summation may be successfully applied to series ()22|l . If we know 
well-defined solutions (like (|28|l is) then Pade summation results may be compared 
with them - it appears that the agreement is unexpectedly good and we expect that 
this occurs for another cases as well. 

A.G. wishes to thank the Polish Ministry of Scientific Research and Information Tech- 
nology for support under Grant No: 1P03B 060 27. 
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